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ABSTRACT 

We analyse CP-violating effects in Z ^ 4 jet decays, assuming the presence 
of CP-violating effective ZbbG and ZbbGG couplings. We discuss the influ- 
ence of these couplings on the decay width. Furthermore, we propose various 
strategies of a direct search for such CP-violating couplings by using differ- 
ent CP-odd observables. The present data of LEP 1 should give significant 
information on the couplings. 
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1 Introduction 



In electron-positron collider experiments at LEP and SLC, a large number of Z bosons has been 
collected so that the detailed study of the decays of the Z boson has been made possible [0]. 
An interesting topic is the test of CP symmetry in such Z decays. There is already a number 
of theoretical and references therein) and experimental ^y]-|2^ studies of this subject. 



In the present paper we will study a flavour-diagonal Z decay where CP-violating effects within 
the Standard Model (SM) are estimated to be very small 0]. Thus, looking for CP violation 
in such Z decays means looking for new physics beyond the SM. 

For a model-independent systematic analysis of CP violation in Z decays we use the effective 
Lagrangian approach as described in [^^. Of particular interest are Z decays involving heavy 
leptons or quarks. Thus, the process Z bbG which is sensitive to effective CP-violating 



couplings in the ZbbG vertex has been analysed theoretically in [^,0 and experimentally in 
12^ . No significant deviation from the SM has been found. 

Here we present an analysis of the 4 jet decays of the Z boson involving b quarks. If CP- 
violating couplings are introduced in the ZbbG vertex, they will, because of gauge invariance 
of QCD, appear in the ZbbGG vertex as well. But the ZbbGG vertex could in principle contain 
new coupling parameters. The 4 jet analysis looks into both, 4- and 5-point vertices. 

In this paper we present the results of our calculations of the process Z ^ 4 jets including 
CP-violating couplings, with at least two of the jets originating from a 6 or 6 quark. The 
following three subprocesses contribute to the 4 jet decay: 

(1) 
(2) 
(3) 

(g = u, d, s, c) . 

We will always assume unpolarized e"^, e~ beams and show the results for each process 
individually as well as the results for the sum of them. In the experiments, of course, only the 
sum of the three processes can be observed easily. 

In chapter 2 we explain the theoretical framework of our computations. Next, in chapter 3, 
we analyse the anomalous couplings for partons in the final state. First, we discuss anomalous 
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contributions to the decay width. Then, we define different CP-odd tensor and vector observ- 
ables and calculate their sensitivities to anomalous couplings. In order to find out how "good" 
for the measurement of the new couplings our observables are, we compare them to the optimal 
observables. In chapter 4 we study decay width, tensor, vector and optimal observables in four 
different scenarios for an experimental analysis. Finally, we compare our results with results of 
the 3 jet decay. Our conclusions can be found in chapter 5. 

2 Effective Lagrangian Approach 

For a model independent study of CP violation in 4 jet decays of the Z boson we use the 
effective Lagrangian approach as explained in |Q. We add to the SM Lagrangian Csm a CP- 
violating term Ccp containing all CP-odd local operators with a mass dimension d <Q (after 
electroweak symmetry breaking) that can be constructed with SM fields. The effective CP- 
violating Lagrangian relevant to our analysis is: 

- '-d, b{x) a^" 75 b{x) [d^ Z,{x) - Z^{x)] 

- '-d',b{x)T^a^''^,b{x)G%{x) 

+ [ hvt b{x) T'^ Y b{x) + hAt b{x) Y 75 b{x) ] Z^{x) , (4) 

where b[x) denotes the b quark field, Z^[x) and G'J^^(x) represent the field of the Z boson and 
the field strength tensor of the gluon, respectively, and T°- = A"/2 are the generators of SU (3)c 
28| . In (^) db is the weak dipole moment, d'j^ the chromoelectric dipole moment, and hyb, 



are CP-violating vector and axial vector chirality conserving coupling constants. As effective 
coupling constants in Ccp the parameters d^, (ij,, hyb, hAb are real. They are related to form 
factors of vertices but should not be confused with the latter (cf. e. g. p!8[]). 



Information on the spin of the final state partons in (|I| - ^ is hardly available experimentally. 
Thus, we consider as observables only the parton's energies and momenta. Then, effects linear 
in the dipole form factors d^ and d'^ are suppressed by powers of rrib/mz- So angular correlations 
of the jets in Z — ^ 4 jets are only sensitive to the couplings hvb and hAb- 

The corresponding vertices following from Ccp are shown in figure |l|. Because the non- 
abelian field strength tensor has a term quadratic in the gluon fields the ZbbG- and ZbbGG- 
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Figure 1: The CP-violating vertices. 



vertices are related. 

We define dimensionless coupling constants hvb,Ab using the Z mass as the scale parameter 

by 

hvb,Ab = -. Z 2" ^Vb,Ab ■ (5) 

sm vw cos Vw rriz 

For numerical calculations we set uiz = 91.187 GeV, sin^ t^vf = 0.2236 and the fine structure 
constant at the Z mass to a = 1/128.9 Our calculations are carried out in leading order 



of the CP- violating couplings of Ccp and the SM couplings. A non- vanishing b quark mass of 
4.5 GeV is included ; masses of u, d, s, c quarks are neglected. 

3 Study of CP- violating couplings for partons in the final 
state 

In this chapter we discuss an ideal experiment where one is able to flavour-tag the partons and 
measure their momenta. We present a study of our CP- violating couplings for each process (|1|) 
- separately and for the sum of them. We have computed the differential and integrated 
decay rates using FORM [^| and M [^| for the analytic and VEGAS for the numerical 



"^We use here the pole mass vahie for the b quark. In our leading order calculation we could as well use the 
b mass at mz'. mi,{mz) — 2.83 GeV [pO|. This would result only in minimal changes in our correlations. 
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calculation. We write the squared matrix element for each subprocess with final state p = 
bbGG, bbbb, bbuu, bbdd, bbss, bbcc in the form: 

+ (hvb + hAb)Sf\<P) + (hvb - hAb)Sf\<P) + hvbhAbSf\ct>) . (6) 



Here stands collectively for the phase space variables, Sq denotes the SM part and 

hb = hAbdYb — hvbdAb , (7) 

hb = hvbgvb — hAbQAb , (8) 

1 2 2 n 1 
gvb = -- + -sm i9w , gAb = --■ (9) 

In the following we drop the index p if the given formula holds for the subprocesses and for 
the sum of the subprocesses. 

The results within the SM have been compared analytically to calculations for vanishing b 
quark mass [^J35|l and to calculations for non- vanishing b quark mass [^. Our results agree 



with these calculations. 

The definition of a 4 jet sample requires the introduction of resolution cuts. We use JADE 



cuts requiring 



with the angle between the momentum directions of any two partons {i ^ j) and Ei, Ej 
their energies in the Z rest system. The expectation value of an observable 0{(f)) is then defined 

as 

jR{<P)d(t) ■ 

3.1 Anomalous contributions to the decay widths 

The solid curves in figure ^ show the results of our calculations for the SM decay widths F'^*^ as 
function of the jet resolution parameter Ucut for the different processes. To check our calculations 
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we computed T^^^ also with the program COMPHEP |38| and found — within numerical errors 
— complete agreement. 
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Figure 2: The decay width for different subprocesses as function of the jet resolution 
parameter y^ut The results for Z bbcc (bbss) are the same as for Z bbuu 

(bbdd). 
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Figure 3: The different contributions to the 4 j^t decay width as function of the jet 
resolution parameter for the sum of the processes (|7| - |^. 



As the decay width is a CP-even observable the contribution of the CP- violating interaction 



to it adds incoherently to the SM one |jT5[: 



r = r^^'^ + AF^^ 



(12) 



with AF being quadratic in the new couplings. In figure g the dashed curves represent AF 
as function of ycut assuming hyt = ^Afe = 1/ v^- 
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As we can see, the dominant decay is (|l|). In comparison to this process, the processes (|^) 
give only contributions at the per cent level, process (H) at the per mille level to the decay 
width. From @ we find: 

Ar^^ = {hvb^ + hAb^)T^ + {hvb - hAb )T^ . (13) 

Because in (^) turns out to be odd under the exchange of quark and anti-quark momenta, 
its integral vanishes. In Figure we compare for the sum of the processes ([l| - 15) F*^*^, F3 and 
F4. For hvb, hAb of order one AF'^'^ is only a correction of a few per cent to F"^*^ even if all 
processes (|l| - 1^) are added up. Thus, considering the theoretical uncertainties in the SM 4 jet 
decay rate, a determination of the new couplings by measuring the decay width alone does not 
look promising. 

3.2 CP-odd observables 

3.2.1 Tensor and vector observables 

We now turn to a study of our CP-violating couplings using CP-odd observables constructed 
from the momentum directions of the h and h quarks, k;, = kfe/|kb| and kj; = kg/lkjl (cf. 

T^f = (ks - k,), (ks X kfe), |k-, X k,r-2 + (, ^ , (14) 

y/") = (k-, X k,), |k^ X kfcp-^ , (15) 

with i, j the Cartesian vector indices in the Z rest system and n = 1,2, 3. 

The observables T^^^ transform as tensors, V^^^ as vectors. For unpolarized e"'"e~ beams 
and our rotationally invariant cuts (|l^) their expectation values are then proportional to the Z 
tensor polarization 5*^ and vector polarization Si, respectively. Defining the positive 2;-axis in 
the e"*" beam direction, we have 





(16) 



(17) 
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where 

53 = '^aVeQAe ^ q , (18) 

9ve + 9Ae 

with gve = — 1/2 + 2sm'^{}w and g^e = — 1/2 the weak vector and axial vector Zee couphngs. 
This shows that the components T33 ■* and V^^^ are the most sensitive ones. 

Note that the tensor observables do not change their sign upon charge misidentification 
(k^ k;,) whereas the vector observables do. Thus, it is only for the measurement of the latter 
that charge identification is indispensable, which makes the vector observables less valuable for 
the experimental analysis. 

We have computed the expectation values of the observables (0), (0) for different JADE 
cuts (p!OD, as function of hb (^ and hb (H). The expectation value of a CP-odd observable O 
has the following general form: 

<0>=(ci/i, + C2M-^, (19) 

J- 4 jets 

where I'fYets T^j^ts denote the corresponding Z — 4 jets decay widths in the SM and in the 
theory with SM plus CP-violating couplings, respectively. In an experimental analysis Tf^ets 
should be taken from the theoretical calculation, r4jets and < O > from the experimental 
measurement. The quantity <0> ■ T4^jets is then an observable strictly linear in the anomalous 
couplings. 

From the measurement of a single observable (|T^ we can get a simple estimate of its 
sensitivity to by assuming hb = 0. The error on a measurement of hb is then to leading order 
in the anomalous couplings: 



^hb = r- , (20) 

Ci V iV 



where N is the number of events within cuts. Similarly, assuming /if, = we get the error on 
hb as 

= I I • ^21) 

A measure for the sensitivity of O to hb (hb) is then l/6hb {l/6hb). However, since we want to 
estimate 2 anomalous couplings /if,, hb we should consider 2 linearly independent observables 
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Oi^2 such that: 

YSM 

J- 4 jets 
YSM 

<02>={hC2l + hC22)^^. (22) 

J- 4 jets 

The sensitivity of these observables to the anomalous couphngs is estimated in the standard 
way. Neglecting terms quadratic in the anomalous couplings the combined measurement of 
< Oi> and <02> with a data sample of N events (within the considered cuts) leads to an 
error ellipse 

{6hbf V{h)^^ + 2 6hb6hb V{h)^^ + {6hbf V{h)^i = 1 . (23) 

Here V(/i) denotes the covariance matrix of the estimated couplings. We have in matrix nota- 
tion: 

y{h) = ^c-'y{o){c~Y , (24) 

cii cu \ ^25) 

^ C21 022/ 

where 

= Js^^ 1 0.{<p)0,{<p)So{<P)d<p (26) 

are the elements of the covariance matrix of the observables Oi, calculated in the SM. A 
measurement of hb, hb has to produce a mean value point outside the ellipse ( pSj) to be able to 
claim a non-zero effect at the 1 s. d. level. 

3.2.2 Optimal observables 

In addition to the tensor and vector observables (0, ^) we study optimal observables, which 
have the largest possible statistical signal-to- noise ratio ||39| , |40|j4l[| . Neglecting higher orders in 
the anomalous couplings the optimal observables for measuring hb and hb are obtained from 
the differential cross section (||) as 

0.(0) = !^, (^ = 1,2). (27) 
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The expectation values for the optimal observables are as in (|2^) with the coefficient matrix 
elements 



1 /-S-mW.,. (28) 



For optimal observables we have 



/ Sodcj) J So{4>) 



= ViO),, , (29) 



Vih)T^' = Nc, . (30) 



3.2.3 Numerical results 



We have calculated the sensitivities to ht and hb for different tensor, vector and the optimal 
observables varying the jet resolution parameter i/cut ■ We assume a total number of 10^ 4 jet 
events from (|l]- ^) for i/cut = 0.01: 

Niy,^t = 0.01) = 10000 . (31) 

The number of events for other values of ycut and for the various subprocesses is then calculated 
within the SM. The total number of Z decays corresponding to (|31|) is Nfot — 6.4 ■ 10^. 



For the process Z bbGG, we found that, in very good approximation, the tensor ob- 
servables are only sensitive to hf, and the vector observables only to hf,. The sensitivities to 
these CP-odd couplings as calculated from ( pOD and (0), respectively, are shown in figure ^ 



The sensitivity decreases with increasing ycut for all observables due to the decrease in number 
of events available. The differences due to the different weight factors for tensor and vector 
observables Tga'', V^''"^ {n = 1,2,3) are only small and all observables considered have nearly 
optimal sensitivities. 

For the processes Z bbbb, Z bbuu and Z bbdd we present plots analogous to those 
for Z bbGG in figures ^ - However, here the correlation of the sensitivity of hh (hb) with 
tensor (vector) observables no longer holds. In Z —>■ bbdd, for instance, the vector observables 
are more sensitive to hf, than to hf,. In the bottom of figure |^ we get a singularity for 6hb 
because at ycut ~ 0.06 the expectation values of our vector observables become zero. Left from 
the singularity all expectation values are negative, right from the singularity they are positive. 
For the sum of the processes the singularity in hf, vanishes, because we have to add up both the 
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variances of all subprocesses in the denominator and the expectation values of all subprocesses 
in the numerator of 

Two results are striking: 

1. The sensitivities obtainable with optimal observables from the subprocesses (0, ^ are 
as good as or even better than those from (P even if the number of events from |^) 
represents only a small fraction of those from (|I]) (cf. figure |^). 

2. For the processes (0, ^ the tensor and vector observables do not reach optimal sensitivi- 
ties. 

One may understand these points in the following way: In the Feynman graph giving the 
CP-violating amplitude for process (^), the gluon which comes out of the CP- violating vertex 
splits into a, q q pair. This means that the CP-violating vertex can be analysed not only by 
using the angular correlations of h and h quark, but also, by means of the momentum directions 
of q and q quarks. If the tensor and vector observables ([l^, p!5|) are used for the measurement 
all information on the CP-violating vertex delivered from the angular correlations of q and q 
quarks is lost but, on the other hand, it is retained by the optimal observables. For process (|^) 
this argumentation is similar. 

In figure ||, results for the sum of the subprocesses (|l| - ||) are shown for ycut = 0.02 
for a combined measurement of T33 and (|1^) and for the optimal observables. The 

comparison of the solid bands (measurement of T^l^ or V^'^'^ alone, respectively) with the dashed 

~ ('2) 
lines corresponding to = ±1 and hi, = ±1 shows that T33 is mostly sensitive to hb and that 

(2) ~ 
V3 is mostly sensitive to h^. The comparison of the inner- and outermost ellipses shows that 

tensor and vector observables do not reach the optimal sensitivity for the sum of the processes (|I] 

-^). This is remarkable since for the dominant process (|I|) they do. Thus, as already discussed 

above, (|^) and (^) which contribute little in the decay rate have a much larger influence in 

CP-odd observables. In tables |l|, |^, |^ of appendix |^, we list the elements of the coefficient 

matrix (^Sj) and the covariance matrix (H) of the observables Oi = T^f (0), O2 = V^^ (||) 

and the coefficient matrix elements (|28|) for the optimal observables for different values of 

the jet resolution parameter i/cut (0)- 
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Z ^ bbGG 




Figure 4: The inverse sensitivities of tensor, vector and optimal observables to hi, 
and hh obtainable in Z —>■ bbGG, as function of the jet resolution parameter i/cut 

^n\ ) assuming (j^Tjj for the number of events. 
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Figure 5: The inverse sensitivities of tensor, vector and optimal observables to hi, 
and hh obtainable in Z —>■ bbbb, as function of the jet resolution parameter y^ut 

( |7Z| j assuming (j^Tjj for the number of events. 
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Z bbuu 




Figure 6: The inverse sensitivities of tensor, vector and optimal observables to hi, 
and hh obtainable in Z ^ bbuu, as function of the jet resolution parameter i/cut 

( |7Z| j assuming (j^Tjj for the number of events. The results for Z —>■ bbcc are identical. 
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Z ^ bbdd 




Figure 7: The inverse sensitivities of tensor, vector and optimal observables to hi, 
and hh obtainable in Z ^ bbdd, as function of the jet resolution parameter ycut 

^n\ ) assuming (j^Tjj for the number of events. The results for Z — * bbss are identical. 
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Figure 8: Contour plot for i/cut = 0.02 of the 1 s. d. errors on hyt, hAh as they 
can be obtained ideally from the measurement of different observables for Z —* 4 j^^s. 
The innermost ellipsis is the result for the optimal observables Oi and O2 at parton 
level. The outermost ellipsis is the result for the combined measurement ofT^, Vg*-^^ 
at parton level. The solid bands correspond to a measurement of T^ or V^"^^ alone, 
respectively. The overlap region of the tensor and vector errors is marked grey. The 
dashed lines correspond to hb = ±1 and hi, = ±1 as indicated. The ellipsis in the 
middle is the combined result for the optimal observables Oi and O2 for analysis 1 
(cf. chapter^. 
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4 CP-violating observables for jets 



In an experimental analysis one can only measure jets as the "footprints" of the underlying 
partons, but not the partons themselves. So if we want to compare our calculations directly to 
experimental data, we must define observables for jets. In LEP experiments it is possible to 



tag a jet referring to a quark with b flavour [53]. In principle one can even distinguish between 



b and b by measuring the jet charge, but this is difficult in practice. We propose four different 
types of analyses with the 4 jet data sample: 

• Analysis 1: One jet comes from b fragmentation, another from b fragmentation (double 
b tag); the other two jets (jets 3 and 4) are ordered according to the magnitude of their 
momenta. 

For the next three analyses, we propose to make an ordering of all four jets according to the 
magnitude of their momenta: 

|qi| > |q2| > Iqsl > |q4| • (32) 

In the following we call jet 1 the jet with the highest magnitude of momentum, jet 2 the jet 
with the second highest magnitude of momentum and so on. 

• Analysis 2: Jet 1 comes from 6 or 6 fragmentation. 

• Analysis 3: Jet 2 comes from 6 or 6 fragmentation. 

• Analysis 4: No requirement to the jet flavour (flavour blind case). 

In appendix ^ we list the different classes of events for each of the subprocesses (|| - 1^) as they 
contribute to these analyses. 

In analyses 2 - 4 we do not distinguish between b and b jets. It turns out that this in essence 
eliminates the dependence of the distributions on the CP-odd parameter hi,. Thus here we can 
only measure hb and we set hb = for these analyses. 

4.1 Anomalous contributions to the decay width 

We computed the total decay width for the 4 jet decays of the Z boson with at least two jets 
coming from b or b fragmentation for the different analyses. Because a momentum ordering of 
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jets can't influence a decay rate, the analyses 1, 4 give results identical to those for the partons 
in the flnal state. In analyses 2, 3 some events are rejected as can be seen from the tables ^ 
10| , \n\ of appendix p. The decay width must decrease in comparison to the other two analyses. 
Figure ^ shows this effect both for the SM contribution and for the contribution of the CP 
violating interaction to the decay width assuming hvb = hAb = 
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Figure 9: The 4 j^t decay width as function of the jet resolution parameter for 
different analyses. The three upper curves are the SM contributions, the others are 
the contributions from the CP-violating interaction if we set hvb = hAb = 1/ v^- 
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4.2 CP-odd observables 

4.2.1 Tensor and vector observables 



We found in chapter ^ that the observables T^f^ (|T^) and ^3*'^'' (|l^) were the most sensitive 
ones. The same was found for the 3 jet decays (cf. Jl^)- Thus, from now on we concentrate 
on this type of observables. 

Analysis 1: The tensor and vector observables in this analysis are the same as for partons: 
T^l^ (|T^ and V^"^^ (JT^). All results are identical to the parton case summed over the subprocesses 
(2 - H) of chapter |^. Thus the sensitivity of a measurement of Tgg'* and V^'^^ to hvb,Ab for 
Vcut = 0.02 is obtained from figure |^. For the measurement of the tensor observable, which is 
C-even, we do not need to distinguish between b and b quark. A sufficient selection criterion 
is then that we demand two jets coming from b gib fragmentation. For the measurement of 
the vector observable, which is C-odd, we need to distinguish between jets coming from b or b 
fragmentation. This can be done experimentally by measuring the jet charge. 

Analysis 2, 3, 4: As tensor observable we chose now 

Tf = (qi - qs). (qi x qa), + ^ j) , (33) 

where q^ = qi/|qi|. We computed the expectation values, variances etc. of the most sensitive 
component z = j = 3 of these observables. All results are shown in figure In table | in 
appendix ^ we list the coefficients of the expectation values (|19D for O = T'fr^ (|33|) for different 
values of the jet resolution parameter i/cut for analysis 3. 

4.2.2 Optimal observables 

The optimal observables are given in (127|) , where (p stands for the relevant phase space variables. 
Note that in calculating Sj{(f)) (j = 0, 1, 2) from (|D we have to sum over the subprocesses (|l| - 
^ taking into account how they contribute to the various analyses (tables |^ - pT]) . 

In figure H we show the results for analysis 1 for i/cut = 0.02 in the /lyb-Zi^b-plane. Compared 
to the tensor and vector observables T33 , combined the optimal observables give only a 
marginal improvement now. This is in contrast to the partonic case and shows again that 
a lot of information about the CP-violating couplings is contained in the distribution of the 
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Figure 10: The error (inverse sensitivity) 6hb obtainable from the tensor observable 
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33 I K 41/ from analysis 1 and T'33 



p^ ) from analyses 2 - 4 function of the jet 



resolution parameter ycut iWW assuming (WW for the number of events. 



secondary quark and anti quark in the subprocesses (0, This information is washed out by 
assuming only knowledge of the momentum ordering of the two corresponding jets. We give 
the numerical values for the elements of the coefficient matrix ( pSD for the optimal observables 
(p7|) for different values of the jet resolution parameter y^ut (|10|) for analysis 1 in table |] of 
appendix ^ 

In figure [ll] we show the inverse sensitivities 6hb for the optimal observable Oi (cf. (^71) ) 
in the analyses 2 - 4, as function of the jet resolution parameter. It is interesting to note that 
using the tensor observable ( ^3]) analysis 3 is superior to 2 whereas with optimal observables 
the reverse is true. In table ^ of appendix 0, we list the coefficients of the expectation values 



(pj|) for O = Oi (|27|) for different values of the jet resolution parameter ycut ( P^ for analysis 3. 
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Figure 11: The error 6hb obtainable from the optimal observable Oi from ana- 
lyses i - as function of the jet resolution parameter ycut l \iW assuming i^3l\ ) for the 
number of events. 



4.3 Comparison with the decay Z ^ 3 jets 

Since hb is in essence only measurable with b and b distinction we concentrate on hf, in the 
following as measured with the tensor observables and the optimal observable Oi in analyses 1 
- 4. To compare the sensitivities of these analyses to those from the 3 jet analyses we calculate 
for each observable O the total number of Z events needed to measure hb with a 1 s. d. accuracy 
6hb within the cuts considered. In figures [1^, |1^ we show these results for analysis 1 and 3, 
respectively. Our results for the 3 jet analyses agree with the calculations |T3|,[T^. We see that 



the 4 jet analyses are competitive and even better than the 3 jet analyses for small values of 
the cut parameter ycut- It should be noted, however, that our results concern the statistical 
errors only. Taking into account experimental efficiencies and systematic errors could change 
the situation considerably. 
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Figure 12: Comparison of the sensitivities of the best tensor observable ( |7^ and 



for the Z ^ 4 J^^s ( analysis 1 ) with the best tensor 
• 3 jets analysis fcf.^F^]). The results for \hb\'^Ntot are 
shown as function of the jet resolution parameter. Ntot is the total number of Z decays 
required to see an effect at the 1 s. d. level for given \hi,\. 



the optimal observable 0\ 
observable i ff^ j for the Z 



5 Conclusions 

In this paper, we have presented various calculations concerning the search for CP violation in 
the 4 jet decays of the Z boson with at least two of the jets originating from h and h quarks. We 
have studied a CP-violating contact interaction with a vector and axial vector coupling hyb, 
hAb (§)• Such couplings can arise at one loop level in multi-Higgs extensions of the Standard 



Model igig. 

We found that, for reasonable values of the coupling constants, the additional contribution 
of the contact interaction to the decay width is at most at the percent level. The decay width 
alone is therefore not appropriate for determining the coupling constants. 

We investigated tensor and vector as well as optimal observables which can be used for the 
measurement of the anomalous couplings. We studied different scenarios for an experimental 
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Figure 13: Comparison of the sensitivities of the best tensor observable T' 



(2) 
33 



and the optimal observable Oi Ml) for the Z ^ 4 j^ts (analysis 3) with the tensor 



observable T''^^^ ^B^) for the Z — » 3 jets analysis (s. also / |7^J7^/ / The results \hb\'^ Nf 



tot 



are shown as function of the jet resolution parameter. Nt^t is the total number of Z 
decays required to see an effect at the 1 s. d. level for given \hh\. 



analysis of the anomalous couplings: The ideal case where all the momenta and flavours of the 
partons can be reconstructed from the jets and four realistic cases where flavour information is 
available only for the h jets. 

If flavour tagging of all jets is available then, with a total number of Z decays Ntot — 6.4- 10^ 
and choosing a cut parameter pcut = 0.02 the anomalous coupling constants hh, hb (0, |) can 
be determined with an accuracy of order 0.1 - 0.2 at 1 s. d. level using optimal observables 
(see figs. I -I). 

In the more realistic case where flavour tagging is available only for b and b jets, the coupling 
constant hb can be measured with an accuracy of order 0.5 - 0.6 using the same total number 
of Z decays. In such a measurement b — b distinction is not necessary. Using in particular the 
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simple tensor observable Tgg"* (14) for the measurement, an almost optimal sensitivity to can 
be attained. 

If b — b distinction is experimentally realizable, the coupling constant hfy can be measured 
with an accuracy of order 0.8 . Again we found a simple vector observable V^"^^ (15) with an 



almost optimal sensitivity to hb- Ub — b distinction is experimentally not realizable the coupling 
constant hf, remains essentially unconstrained from measurements of CP-odd observables. It 
can be bounded indirectly by assuming, for instance, that its contribution to the 4 jet width 
does not exceed 5%. This implies then \hh\^0{l). 

In our theoretical investigations we assumed always 100% efficiencies and considered the 
statistical errors only. But the total number of Z decays collected by the LEP and SLC exper- 
iments together is of order 10^. Thus the accuracies in the determinations of h^, hb discussed 
above should indeed be within experimental reach. 

Comparing 3 and 4 jet analyses we found that the sensitivity to the anomalous coupling hb 
was roughly constant as function of the cut parameter i/cut for i/cut < 0.1 in the 3 jet case. For 
the 4 jet case the sensitivity was found to increase as i/cut decreases. For T/cwi ^0.01 the 4 jet 
sensitivity was found to exceed that from 3 jets (figures 0, 0). Of course in an experimental 
analysis one should try to make both 3 and 4 jet analyses in order to extract the maximal 
possible information from the data. 

For the experimental analyses, one usually has to make Monte Carlo simulations. For this 
purpose one needs matrix elements including the CP-violating interaction. These are available 
from us in the form of FORTRAN subroutines.^ 

To conclude: we have discussed in detail various possibilities to determine or obtain limits 
on anomalous CP-violating ZbbG and ZbbGG couplings. As shown in p6| , ^ this will give 



valuable information on the scalar sector in multi-Higgs extensions of the Standard Model. 
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A Numerical Values 



We list some numerical results for the coefficient matrices and covariance matrices in different 
studies. The statistical errors of the numerical calculation are typically at the per cent level. 



Ucut 


Cll 


Cl2 


C21 


C22 


0.01 


-0.01827 


-2.496 ■ 10-5 


-3.556-10-^ 


6.426 ■ 10-3 


0.02 


-0.02320 


-6.706 ■ 10-5 


-4.504-10-^ 


8.336 ■ 10-3 


0.05 


-0.03253 


-1.666 ■ 10-^ 


-7.127-10-^ 


0.01248 



Table 1: The numerical values for the elements of the coefficient matrix of the 
observables Oi = T33 ([T^j, O2 = ( |73t ) for different values of the jet resolution 



parameter ycut (H^ for partons in the final state (section 



Vcut 


ViO)u 


V{0)i2 = ViO)21 




0.01 


0.2772 


0.01811 


0.1109 


0.02 


0.3429 


0.02334 


0.1427 


0.05 


0.4332 


0.02949 


0.1959 



Table 2: The numerical values for the elements of the covariance matrix ( ^q ) of the 



observables Oi = T33 ({T^), O2 = ( \r^ ) for different values of the jet resolution 
parameter Ucut (IHj for partons in the final state (section [^. 



Vcut 


Cll 


C12 — C21 


C22 


0.01 


5.556 ■ 10-3 


3.228 • 10-3 


4.726 - 10-3 


0.02 


9.557-10-3 


6.698-10-3 


9.588 - 10-3 


0.05 


0.02274 


0.02002 


0.02816 



Table 3: The numerical values for the elements of the coefficient matrix ([^^ for the 
optimal observables for different values of the jet resolution parameter Ucut ( f^j 
for partons in the final state (section 
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Vcut 


Cl 


0.01 


-0.01020 


0.02 


-0.01256 


0.05 


-0.01502 



Table 4: The numerical values the coefficients of the expectation values ( |7^ j for 
O = for different values of the jet resolution parameter i/cut ( T^j for analysis 

3 of section 0. 



Vcut 


Cll 


C12 = C21 


C22 


0.01 


1.514- 10-3 


1.041 • 10-^ 


4.825 • 10-^ 


0.02 


1.981 • 10-3 


1.651 ■ 10"^ 


6.396 ■ 10"^ 


0.05 


2.956 • 10-3 


3.042 ■ 10"^ 


9.907-10-^ 



Table 5: The numerical values for the elements of the coefficient matrix for the 
optimal observables for different values of the jet resolution parameter ycut 
for analysis 1 of section 0. 



Vcut 


Cl 


0.01 


5.535-10-^ 


0.02 


7.124-10-^ 


0.05 


8.627-10-^ 



Table 6: The numerical values for the coefficients of the expectation values ( fT^j for 
the optimal observable O = Oi ^Ij ) for different values of the jet resolution parameter 
Vcut ( P^ j for analysis 3 of section 0. 



B Eventclasses 

Here we explain which classes of events contribute to the four different analyses as defined in 
chapter First, we compare the partonic phase space with the jet phase space. 
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Process 


Phase space restriction 


e+e- -^Z ^ bbGG 


|ki| > Ikal 


e^e~ Z ^ bbbb 


|k+| > |q+|, |k_| > |q_| 


e^e~ Z ^ bbqq 





Table 7: The restrictions on the partonic phase space for the different processes due 
to identical particles in the final state. 



Analysis 


Phase 


space restriction 


1 






> 


q4 




2, 3,4 


|qil > 


q2 


> 


qs 


> |q4| 



Table 8: The restrictions on the jet phase space in the analyses 1-4- 

In Tables 0, | we list the restrictions on the phase space for the partonic processes (|1] - 1^) 
and for the jets in the analyses 1 - 4 as defined in chapter ^. 

In tables ^- |n| we list all possibilities how the 4 partons of the reactions (|T]- §|) can give 4 
jets with the ordering criteria of the analyses 1-4. The full points indicate that an event class 
satisfies the respective selection criterion. 
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6 jet 


6 jet 


jet 3 


jet 4 


Analysis 1 








G{k2) 


• 



jet 1 


jet 2 


jet 3 


jet 4 


Analysis 2 


Analysis 3 


A 1 • A 

Analysis 4 




b{k-) 


G{ki) 


G{k2) 


• 


• 


• 




Gik,) 


6(A;_) 


G{k2) 


• 




• 


KK) 


G{ki) 


^(^2) 


b{k^) 


• 




• 




KK) 


G{ki) 


G{k2) 


• 


• 


• 


b{k-) 


G{ki) 


Kk+) 


G{k2) 


• 




• 


Kk-) 


G{ki) 


G{k2) 


Kk+) 








G{ki) 


Kk+) 


b{k-) 


G{k2) 




• 




Giki) 


b{k+) 


G{k2) 


Kk-) 




• 




G{ki) 


b{k_) 


b{k^) 


G{k2) 




• 




G{k,) 


6(fc_) 


G{k2) 


b{k^) 




• 




G{k,) 


G{k2) 


b{k^) 


b{k_) 








G{k,) 


G{k2) 


b{kJ) 


Kk+) 









Table 9: The only possibility for the partons in e'^e~ — > Z — > bbGG to fulfill the 
selection criterion of analysis 1 and the 12 possibilities for them to give 4 momentum 
ordered jets. It is indicated by dots which event class contributes to which of the 
analyses. 
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6 jet 


6 jet 


jet 3 


jet 4 


Analysis 1 






Kk+) 


b{k^) 


Kq-) 


b{q+) 


• 








b{k-) 




b{q-) 


• 








b{k-) 


b{q-) 


Kk+) 


• 






Ku) 


b{k_) 


b{k+) 


b{q-) 


• 






KK) 


Kq-) 


b{k^) 


Kq+) 












Kq+) 


b{k-) 










Kq-) 


b{k-) 


b{k+) 












b{k+) 


b{k-) 
















jet 1 


jet 2 


jet 3 


jet 4 


Analysis 2 


Analysis 3 


Analysis 4 




bik.) 


Kq-) 


biq+) 








KK) 


b{k_) 




b{q-) 








KK) 


Kq+) 


b{k_) 


b{q-) 








b{k.) 




b{q-) 


Kq+) 








b{k.) 


KK) 




b{q-) 














b{q+) 









Table 10: The 8 possibilities for the partons in e^e~ Z ^ bbbb to fulfill the selec- 
tion criterion of analysis 1 and the 6 possibilities for them to give 4 momentum ordered 
jets. It is indicated by dots which event class contributes to which of the analyses. 
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6 jet 


6 jet 


jet 3 


jet 4 


Analysis 1 




Kk~) 






• 




bik-) 






• 



jet 1 


jot 2 


jet 3 


jot 1 


Analysis 2 


A 1 • 

Analysis 3 


Analysis 4 


Kk+) 


6(A:_) 


QiQ-) 




• 


• 


• 


Kk+) 




q{q+) 


q{q-) 


• 


• 


• 


Kk+) 




b{k^) 




• 




• 


Kk+) 


9(9+) 


b{k^) 




• 
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b{k+) 






b{k-) 


• 




• 


b{k+) 






b(k-) 


• 




• 


b(k-) 
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• 


• 


b{k_) 


6(A;+) 


^(?+) 




• 


• 


• 


h{k_) 




Kk+) 




• 




• 


h{k_) 




Kk+) 




• 




• 


b{k-) 






b{k+) 


• 




• 


b{k-) 






b{k+) 


• 




• 


q{<}-) 






<}(<}+) 




• 


• 






b{k-) 






• 


• 




Kk+) 




6(A;_) 




• 


• 






lii-) 


6(A:_) 




• 


• 




6(A;_) 


Kk+) 


9(9+) 












b{k+) 


q{q-) 










6(A;_) 




b{k^) 










Hk-) 


<li'l-) 












9(9+) 


b{k+) 


6(A;_) 












Kk+) 


6(A:_) 












b{k^) 


Hk+) 












b{k-) 


b{k+) 









Table 11: The 2 possibilities for the partons in e+e~ — > Z — > bbqq to fulfill the se- 
lection criterion of analysis 1 and the 24 possibilities for them to give 4 momentum 
ordered jets. It is indicated by dots which event class contributes to which of the 
analyses. 
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